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Abstract 

Gravitational lensing assists in the detection of quasar hosts by amplifying and 
distorting the host light away from the unresolved quasar core images. We present 
the results of HST observations of 30 quasar hosts at redshifts 1 < z < 4.5. The 
hosts are small in size (r e < 6 kpc), and span a range of morphologies consistent 
with early-types (though smaller in mass) to disky/late-type. The ratio of the black 
hole mass (Mg#, from the virial technique) to the bulge mass (Mb u i ge , from the 
stellar luminosity) at 1 < z < 1.7 is broadly consistent with the local value; while 
Mbh/ Mbuige at z > 1.7 is a factor of 3-6 higher than the local value. But, depending 
on the stellar content the ratio may decline at z > 4 (if E/SO-like), flatten off to 
6-10 times the local value (if Sbc-like), or continue to rise (if Im-like). We infer that 
galaxy bulge masses must have grown by a factor of 3-6 over the redshift range 
3 ^ % ^ 1, and then changed little since z ~ 1. This suggests that the peak epoch 
of galaxy formation for massive galaxies is above z ~ 1. We also estimate the duty 
cycle of luminous AGNs at z > 1 to be ~ 1%, or 10 yrs, with sizable scatter. 
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1 Introduction 



Elsewhere in these proceedings the reader can find summaries of previous work 
on quasar hosts. We concentrate here on the benefits and challenges of using 
gravitational lensing as a technique for measuring the host galaxy, especially 
in regimes where traditional direct imaging methods are difficult - at high red- 
shift, or when the host is either sub-luminous or compact. The context for this 
work is the coevolution of galaxies and supermassive black holes. In the local 
universe, a tight relation is observed between bulge m ass (.Mhnw) and black 



hole mass (.Mrh) measur ed with stellar kinematics ( Gebhardt et al. . 2000; 



Ferrarese &: Merritt . 2000[ ). We now extend the relation to z > 1 using the 
virial technique to estimate M BH (e.g. Kaspi et al., 2000; Vestergaard fc Peter- 



son, 2006), and the stellar lum inosity to estimate .M bulge ( Kormendv fc Richstond . 
19951 Magorrian et al. . 19981 ). The existence, slope, and scatter of a .Mbh/.M bulge 
relation at high redshift can be used to analyze the relative growth rates of 
galaxies and their (presumably ubiquitous) central engines. 

The CfA- Arizona Space Telescope Lens Survey (CASTLES) is a project to im- 
age all known, multiply-imaged, quasars in a homogeneous set of optical and 
near infrared passbands using the Hubble Space Telescope (HST). With the 
number of lens systems now near 100, data are in hand for 80 targets. The ob- 
servations are shallow, 1-2 orbits per filter, but the excellent surface brightness 
sensitivity of the HST leads to a host detect i on in most cases. The overall CAS- 
TLES project is describ ed bv iFa 



host search are given bv lRix et al. 



co et al. (|200lh : early results in the lensed 



( 200l| ) in the same conference proceedings; 



and detect ions and models of the hosts of several individual quasars have b een 
published f)lmoev et all 119981 iKochanek et~ail l2000l : iKeeton et all l2000h . 



Gravitational lensing is a large and growing field of astrophysics so only the 
rudiments can be given here; a number of excellent boo k and reviews are avail- 
able for the full formalism and diverse applications (iBlandford &: Naravan , 
19921: ISchneider. Ehlers. fc FalcoL Il99l IClaeskens fc Surdeil . l2002l) . About one 
in 500 quasars has a sight-line passing close enough to the central potential of 
a massive galaxy for multiple image formation. It has taken surveys of tens of 
thousands of radio and optically selected quasars to yield the sample of ~100 
objects (see the CASTLES web site 1 and that of Liege group 2 ). Even though 
adaptive optics techniques from the ground are improving, stable point spread 
functions, obtained with HST, are still essential for reliable modeling. 



In gravitational lensing, the AGN is magnified into multiple images, but re- 
mains unresolved, whereas the extended light from the host galaxy maps into 
arcs or Einstein rings (ER). A lens model is needed to extract the full infor- 



1 http:/ /www. cfa.harvard.edu/castles 

2 http:/ /vela.astro. ulg.ac.be/themes/extragal/gravlens 
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Fig. 1. An example of the lens modeling technique, (a) Original NIC2 image of 
PG 1115+080 (zq so = 1.72). (b) The host galaxy Einstein ring, after removing the 
best fit lensing galaxy and quasar point sources, (c) The best fitting residuals, (d) 
The parametric model of the host galaxy in the source plane. 



mation content of the lensed host light. In principle, a typical ER having a 
radius ~ 1 arcsec means that HS T imaging potentially yields 50-100 resolution 
elements in the host galaxy in the deep images from the survey. 



2 Modeling 



The image modeling (Fig. 1) uses a custom-built program called LENSFIT 
(Peng et al. 2006, in prep.), which is ba sed on a methodo logy that has been 
well-tested with the GALFIT algorithm ((Peng et al. . 2002j b The model for the 
light profiles of the host and foreground (lens) galaxy uses a Sersic model with 
a concentration index n that is often used to quantify the gross morphology of 
galaxies (e.g. n — 1 for late- type, while n = 4 for early- type). Both the quasar 
point source and the host galaxy light profile are propagated through the lens 
model to produce the image distortion, and multiple images. External shear is 
included to model the tidal influence due to neighbors. All the parameters are 
simultaneously varied to reduce the \ 2 on a pixel-by-pixel basis. The models 
are often very robust in well resolved systems (9 > 1 arcsec), due to the spatial 
separation between the host and the lens, and their different shapes. 
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Fig. 2. The (deprojected) host galaxy luminosity with redshift. Open points: radio 
loud quasar hosts. Solid points: radio quiet or unknown (most likely radio quiet). 
The stellar evolution tracks correspond to zj orm = 5, all normalized to L* by z = 0. 
The E/SO and Sb/c tracks have an initial burst of stars, followed by star formation 
rates that produce colors of E/SO and Sb/c type galaxies by z = 0. The "No 
Evolution" model corresponds to a redshifted spectrum of a z = E/SO galaxy. 



3 Results 



Here we present a summary of our findings, which are detailed along with a de- 
scription of our analysis techniques elsewhere (Peng et al. 2006, in prep.). We 
select on lensing geometry size (6 > 0.7 arcsec), which does not a priori bias 
the intrinsic AGN luminosity selection, or the host luminosities. Therefore, we 
expect the sample of the AGNs to be randomly drawn from the AGN luminos- 
ity function, where the lower limits are determined by various lensing search 
programs. The heterogeneity of different surveys, however, will not affect our 
primary conclusions about the relationships between A^bh an d A^buige, since 
both quantities are measured in the same objects. 



3.1 General Properties 



Figure 2 shows the iJ-band host galaxy luminosities (lensing distortion re- 
moved) versus redshift. Overall, the host luminosities f r om CASTLES appear 
to ag ree well with non- lensing studies ( Kukula et all 12001 : Ridewav et al. 



2001). The host luminosities range from 1 to 20 times L v today, while the 



AGNs are to 3 magnitudes brighter than the host in restframe B to V 
band. Despite their brightnesses, the host galaxies appear to be fairly small in 
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size (typical r e < 6 kpc) for their central AGNs. As we shall see later, when 
coupled with information about their A^bh; both the host luminosities and 
sizes lead to the conclusion that the bulges may be undermassive compared 
to present-day normal galaxies. Lastly, the Sersic index values suggest that 
while a number of quasar hosts at z > 1.5 have steep central concentrations 
consistent with the presence of a bulge, many (30% - 50%) also have low Sersic 
values (n < 2) more analogous to later-type galaxies today. Even those galax- 
ies with high Sersic indices may not qualify as bona fide - fully formed and 
passively evolving - ellipticals, given their small sizes and black hole masses. 

Radio-loud (RLQ) vs. Radio-quiet (RQQ) hosts Quantifying differences be- 
tween RLQ and RQQ hosts has historically been controversial. Figure 2 shows 
that there is not a clear difference in the host luminosities between RLQ (open) 
and RQQ (closed) AGNs in the lensing sample. The diverse selection criteria 
from disjointed surveys are, however , hard to quantify. It is worth to keep 



in mind, however, that in a study bv iKukula et al.l ([20011 ) RLQs were drawn 



from the rare and extreme radio-loud sources which may require atypically 
large BHs to produce. Consequently, the finding of luminous hosts in those 
RLQs may reflect a correlation between A^bh and A^buige at high redshifts. 
The issue of radio correlation with host properties remains unsettled. 

AGN Duty Cycle We can estimate the duty cycle of nuclear activity for 
each object with a host detection. A rough estimate of the duty cycle is: 
D ~ &q(Lqso, z) /^ciLcai, z), where $q and $q are the luminosity functions 
of quasars and galaxies, respectively, appropriate to a given redshift. At z > 1, 
the median duty cycle is 1%, or 10 7 years, with a sizable scatter. 



3.2 Black Hole vs. Bulge Evolution 



Based on quasar and host luminosities we can study the A4bh vs. bulge prop- 
erties at z > 1 (see Peng et al. 2005 for details), where AIbh's are obtaine d 



using the virial technique ((Kaspi et all l2000t IVestergaard. fc Peterson! . |2006) 



1-7 < z < 4.5 Fig. 3 shows the AIbh vs. the restframe i?-band bulge 
luminosity (Lr) for the host galaxies at 1.7 < z < 4.5. Determining Lr 
requires i^-corrections, computed using an Sbc SED; the dashed lines shows 
the small systematic effect of using an E/S0 (right) or Im (left) SED. It is clear 
that a correlation between AIbh and bulge luminosity was already present at 
a lookback time of 10-12 Gyr. Remarkably, the high-redshift hosts appear to 
lie on the same relation as z = normal galaxies, implying that the high- z 
hosts are undermassive in comparison. To explain why, Fig. 3b shows the host 
luminosities after we account for passive evolution of dM^/dz = —0.8 mag; 
specifically, we hold M BH fixed and shift the color points in Fig. 3a to the left. 
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Fig. 3. The relationship of the black hole mass, Mbyl, vs. bulge absolute luminosity 
(Lr, bottom axis; Mr, top axis), at low z (solid round points) and z > 1.7 (open 
points). Solid lines: fit to z « solid points. All open points assume a modern-day 
Sbc-type SED for i^-correction and their average is represented by dotted lines. 
Dashed lines illustrate assumptions of bluer(Im, left) /redde r(E, right) SED s . Ope n 
circle: gravita tionally lensed quasa r hosts. Open triangles: iRidewav et al. 
Open squares: iKukula et al J ([200 ll ). Vertical line in points: a possible lower limit in 
A^bh due to AGNs being broad absorption line QSOs. Criss-crossed points: poten- 
tial problem with lens identification, host detection, radio-loud quasar, or narrow 
line AGN. Panel (a): The observational data. Panel (b): The same data in a, but 
the open points are shifted horizontally by ass uming that the ho sts evolve passively 
with Zf = 5 by dM^/dz = —0.8 mag. See also Peng et all (|2005h for details. 



Now the z > 1.7 hosts are displaced from the local host relation (solid line) by 
a factor of 3-6 in luminosity, which translates into a mass deficit of a factor 
of 3-6 in the quasar hosts compared to local galaxies with the same .Mbh- 



1 < z < 1.7 In contrast, by z ~ 1, Fig. 4 shows that the mass deficit of the 
hosts is mostly reduced (to within a factor of unity in mass), after accounting 
solely for passive evolution. Thus, massive bulges that correspond to luminous 
E/SO galaxies today may have been nearly assembled by z ~ 1. 

By requiring that the hosts evolve onto the local .Mbh vs. Lr of Figs. 3 and 
4, we can illustrate a growth in the M.-qk/ M- bulge ratio with redshift relative 
to today, shown in Figure 5. We find that the ratio of A^bh/A^ bulge increases 
roughly to a factor of 10 higher than today (assuming an Sbc-type SED) out at 
z = 4. This conclusion depends somewhat on the assumption of the SED and 
evolutionary history: the ratio at earlier times would be higher than shown for 
a bluer SED than Sb/c, or for a faster fading rate than passive evolution. 



6 



M R (bulge) [mag] 



11 
10 
51 9 



CD 



6fl 



8 r 



-2 7 



6 r 




-20 



-25 



i i i | i i i i 

(b) z form = 5 + 
pass. evol. 




Oz<1.7 CASTLES 
□ z»l Kukula 
• z«0 Normal 



10 11 12 9 10 11 12 
lo gio( L R/ L o) ( bul S e ) 

Fig. 4. The same diagram as Figure 3, except for redshift of 1 < z < 1.7 quasar 
hosts. The vertical line in the square points indicates that the AGN either has a 
strong narrow Mg II line component or is strongly absorbed in the wings, causing 
a potentially low AIbh estimate. The dotted line is displaced from the solid line, 
representing the local Mbr-Lr relation, by —0.5 (Fig. a) and +0.5 (Fig. b) magni- 
tude. Note the very slight bias between the non-lensed and lensed datasets, which 
might be explained by the difference in the median redshift of (z) me< j = 1-45 for the 
lens sample and (z) me( j = 0.94 for the non-lenses. However, three of the non-lensed 
data points may also have a lower limit on the AIbh estimate, as noted above. 
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Fig. 5. The growth of the A^bh/AI* ratio as a function of (a) redshift and (b) age of 
the universe in Gyrs. Circles are gravitational lens data poin ts, while triang l es are 
from direct imaging of hosts using HST NICMOS #-band feidgwav et all l200ll : 
Kukula et all 1200 if ). Point styles are the same as Fig. 3. The AIbh/AI* ratio ap- 
pears to rise quickly beyond z ~ 1 and may slow, and perhaps flatten, to a factor 
of 6 — 10 x local value by z ps 3. A fading rate of dMu/dz = —0.8 is assumed here 
(passive evolution since zj orm = 5). 
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4 Conclusions 



Detailed modeling of 30 well-observed systems from a total sample of 80 lensed 
quasars has provided new insights into the properties of host galaxies at 1 < 
z < 4.5. About half have Sersic model fits indicative of early type galaxies. 
However, combined with their small sizes of r e < 6 kpc, luminosities, and 
A^bh, it appears that luminous, fully- formed, ellipticals are in a minority 
as hosts of luminous quasars at z > 2. No difference is seen between the 
luminosities of radio-loud and radio-quiet quasars in the sample, with a caveat 
on sample selection. Even at z > 2, the host galaxies follow nearly the same 
relationship between A^bh and luminosity as at low redshifts, but the bulges 
must gain in mass by a factor of 3-6 between 1 < z < 4.5. However, by 
z ~ 1, the mass deficit is mostly gone. Thus massive bulges at z ~ 1 may be 
consistent with being passively evolving, or may still grow by at most a factor 
of 1. Our estimate of the AGN duty cycle is m 1%, or 10 7 years. Ongoing 
work includes using color information to constrain the host star formation 
histories, obtaining sub-millimeter data to measure star formation rate, and 
characterizing detailed host morphology with deeper HST imaging. 
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